Heat capacity of mesoscopically disordered superconductors: implications to MgB 2 



O 

o 

(N 

o 

Q 



A. M. Gabovich 1 , A. I. Voitenko 1 , Mai Suan Li 2 , H. Szymczak 2 , and M. Pekala 3 
1 Crystal Physics Department, Institute of Physics, NUAS, Prospect Nauki 46, 03022 Kiev-22, Ukraine 

2 Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland 
3 Department of Chemistry, University of Warsaw, Al. Zwirki i Wigury 101, PL-02-089 Warsaw, Poland 

Electronic heat capacity C(T) was calculated for a mesoscopically disordered s-wave superconductor 
treated as a spatial ensemble of domains with a continuously varying superconducting properties. 
The domains are assumed to have sizes L > £o, where £o is the coherence length. Each domain 
is characterized by a certain critical temperature T c o in the range [0, T c ]. The averaging over the 
superconducting gap distribution leads to C(T) ~ T 2 for low T, whereas the specific heat anomaly 
at T c is substantially smeared. The results explain well the C(T) data for MgEb, where a multiple- 
gap structure is observed. 
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PACS numbers: 74.20.-z, 74.25.Bt, 74.80.-g 



An un- 
expected discovery of relatively high-temperature, high- 
T, superconductor MgB2 with a critical temperature 
T c « 40KEJ have cast a considerable doubt on the va- 
lidity of the opinion that high T c 's are appropriate to 
substances with a spin-fluctuation-driven Cooper pair- 
ing and, consequently, with a predominantly ci x 2_„2-wave 
symmetry of the superconducting order parameter. In- 
deed, an obvious absence of magnetic ions, a consider- 
able isotopic effecta, and Bardeen-Cooper-Schrieffej'-like 
(BCS-like) coherent peaks in optical conductivity^] and 
spin-lattice relaxationa are indicative of the conventional 
s-wave pairing in MgB2. As for the electron-phonon 
background of superconductivity, it also seems highly 
probable, although a fairly exotic multiple-gap scenario 
is needed to reconcile the available data (see, e.g., Ref. 
|5|). It is very remarkable that the multiple-gap conven- 
tional Cooper pairing is directly found in a number of 
point-contact, tunneling and Raman measurements!! It 
is even more important that the distribution of gaps may 
be rather broacfl and spatially-resolveda, although much 
controversy exists over the number and widths of gaps in 
the electron density of states (DOS). 

In actual truth, the order parameter symmetry for 
MgB2 is not unambiguously determined. The low-T, 
asymptotics of the magnetic field penetration depth A(T) 
was shown by muon spin-rotationo and opticaH measure- 
ments to be a power-law one. This was interpreted as 
either an unconventional superconductivity or at least 
as a highly anisotropic s-wave pairing. Thermodynamic 
measurements might be decisive in determining the low- 
T symmetry-based superconducting properties of MgB 2 
because the minority phases or grain boundaries do not 
affect the results substantially in contrast to, e.g., trans- 
port phenomena. The electronic heat capacity, C(T), 
behavior near T c is also of great importance to elucidate 
the nature of superconductivity here. And, indeed, there 
were a lot of specific heat-.investigations for MgB 2 per- 
formed by various groupsL3. 

The main features of the data for C(T) are (i) small 



valueg-pf— j;he phase transition anomaly 



C s — C n 

at T c ll9llij in comparison to the BCS caseE-i, when the 
ratio n = AC/[-f S {T c )T c } is equal to mbcs = 12/[7C(3)]; 
and (ii) deviations from the asymptotic BCS behavior at 
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Here 7s is a Sommerfeld constant, the subscripts s and 
n correspond to the superconducting and normal states, 
respectively, N(Q) is the electron DOS at the Fermi level, 
Ao is the energy gap value at T = 0, Hb = % = 1. 
The deviations from Eq. (Q) may be twofold • power-law- 
like ~ T 2 B and of the form ~ exp(-^)EE3, where the 
constant A is much less than —T r w 1.76 T r , as it should 



to 



be in the weak-coupling superconductor^, 7 = 1.78. 
is the Euler constant. Thus, the raw specific heat data 
do not give definite answers to the problems of the order 
parameter symmetry and the underlying mechanisms of 
superconductivity. 

In this article on the basis of the experimentally proved 
distribution of energy gaps we show that both main fea- 
tures of C S (T) can be explained by the conventional s- 
wave superconductivity, so that these data can be easily 
reconciled with other observationgJu. The adopted ap- 
proach, being the outgrowth of the earlier onet£l, is phe- 
nomenological because the origin of the gap distribution 
is not knowja precisely. However, in accordance with tun- 
neling dataO, the gap distribution is considered to occur 
in the real space rather than in the k-space, as was sug- 
gested, e.g., in Refs. The theoretical description 
of such spatially disordered superconductors depends on 
the ratio between the characteristic superconducting do- 
main size L and the coherence length £0"- If L > £q, 
superconducting properties are determined by local val- 
ues of the order parameter A. In the case of MgB 2 there 
is a large scatter of £0, inferred from. different experiments 
and for different kinds of samplesO, so that we may es- 
timate this quantity as lying in the range from 25 A to 
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65 A. This dispersion of £o qualitatively correlates with 
the bro&dr-SiPiectra of gaps in tunnel and point-contact 
spectra&ElHa 

Let us examine a T-indepcndcnt configuration of meso- 
scopic domains, with each domain having the following 
properties: 

(A) at T — it is described by a certain superconduct- 
ing order parameter Ao < A™ ax ; 

(B) up to a relevant critical temperature T c0 (A ) = 
■JAq, it behaves like an isotropic BCS superconductor, 
i.e. the superconducting order parameter AXT) is the 
Miihlschlegel function A(T) = A BC s(A , T)H, and the 
electronic specific heat is characterized in this interval by 
the function C S (A,T); 

(C) at T > T c o it transforms into the normal state, 
and the relevant property is C n (T). 

At the same time, the values of Ao scatter for vari- 
ous domains. The current carriers move freely across do- 
mains and inside each domain acquire appropriate prop- 
erties. The adopted picture is especially suitably for su- 
perconductors with small coherence lengths £qi3. 

In other words, each domain above its T c q is in the 
normal phase and its specific heat iai-a 



(C n (T)) = C n (T)p n (T), 



C n (T) 



-N(0)T. 



(2) 



For simplicity we restrict ourselves to the situation when 
the whole sample above T c is electronically homogeneous, 
i.e. is characterized by a common approximately con- 
stant A^(0) value. Below T c0 for a given mesoscopic do- 
main, a corresponding isotropic gap appears on the Fermi 
surface. The microscopic background of the assumed 
scatter in T c o's may be diverse but ultimately manifests 
itself as a variation either of the electron-phonon inter- 
action magnitude or of local values of the Coulomb pseu- 
dopotential. 

In the framework of our phenomenological approach, 
superconductivity (if any) inside a chosen domain is de- 
scribed by the relevant parameters Ao and T c q. They are 
bounded from above by Ao° ax and T Cl respectively. These 
Ao's may or may not group around a certain crowding 
value Ao depending on the sample texture. The existence 
of such two possibilities is in accordance with the varied 
data for MgB2Hi3'E2rBlla. The specific gap distribution is 
described by the function /o(Ao). 

Thus, for all T in the interval [0,T C ], where T c = 
max T c o, the superconducting sample consists of super- 
conducting (s) and nonsuperconducting (n) grains more 
or less homogeneously distributed over the sample vol- 
ume. 

The measured C s (T) is an averaged sum of contribu- 
tions from both phases 



(C(T)) = (C„(T)) + (C S (T)) 



(3) 



which depends on the distribution function /(A, T) of 
superconducting .domains, and on the fraction p n (T) of 
the normal phasdl3 



(C.(T)) = 



, (T) 



C s (A,T)f(A,T)dA. 



(4) 
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Here A max (T) = A B cs (A najc ,T') and /(A, T) is a result 
of the thermal evolution of the initial (at T = 0) distri- 
bution function /o(Ao). It is convenient to normalize all 
temperatures by T c and all energy parameters by Ao™ ax : 
t = T/T c ,5 = A/Ajf ax with relevant indices retained, 
and to consider C S {T) and C n {T) together with their 
averaged counterparts, normalized by the C n (T c ) value, 
i.e., c Stn (t) = C Sin (T)/C n (T c ). Then one can easily find 
that for each domain, characterized by the parameter 6q 
at t = 0, the dimensionless heat capacity is either 



t (t) = t,t>S 



c s (t) = S cbcs ( j- ) , t < S , 



(6) 
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where cbcs {%) is a well-known normalized heat-capacity 
function for a standard BCS superconductorllJ. For a sur- 
mised domain ensemble a distribution function f(A,T) 
for finite T is defined by the formula 

/(A,T)dA = /o(A )dA . (8) 

Then the dimensionless heat capacity takes the form 



(c s (t))=J c BCS i^-j f Q (S a )S dSo. (9) 

Introducing a new variable z = t/Sg and expanding the 
function /o (t/z) into a series we arrive at the proper low- 
t asymptotics 

(c s (t - 0)) = t 2 / -/o (0) c BCS (z) « 2.45 /„ (0) t\ 
Jo z 

(10) 

The i-dependence of the next term in the expansion for 
(c s (t)} can be estimated in the limit t — > by substitution 
of the normalized expression (Q) for cbcs ( z )- It turns 
out that this expression decreases as 0[t 5 ^ 2 exp(— ^)]. 

Now, in the same low-T region let us take a look at 
the contribution (c n (t)) of the continuously expanding 
normal phase. At any T, all domains with Ao < (i.e. 
Sq < t) are nonsuperconducting, with the total normal 
phase fraction being 



Pn{t) = Pn(Q) 



f (S )dS . 



(11) 



For simplicity, below we restrict ourselves to the case 
when all domains at t = are superconducting, i.e. 
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Pn(0) = 0. A generalization to the case p„(0) / 
is obvious: at each temperature there exists an addi- 
tional contribution from the normal phase. Then the 
function fo(So) should be normalized by 1 — p n (0), and 
all averaging-driven effects would accordingly decrease. 
Moreover, if p n (0) ^ 0, the observed heat capacity (c(t)) 
must include an extra linear contribution p n (0)t in the 
true superconducting state exhibiting the Meissner effect. 

As for the second term in Eq. (11), the approximation 
of /o(<5o) by its limiting value /o(0) demonstrates that 
the main temperature-dependent contribution to p n (t) 
is linear in t. Since c n (t) is also a linear function of t 1 
the apparent contribution (c n (t)} of the normal phase 
to the resulting specific heat (c(t)) is quadratic in t for 
small t, similarly to (c s (t)). Thus, in the suggested model 
of the disordered superconductor with a broad continu- 
ous spatial distribution of domains, possessing different 
T c 's, normal and superconducting contributions to ther- 
modynamical quantities are functionally indistinguish- 
able from each other. 

In addition to the low-T asymptotics the overall In- 
dependence of the heat capacity C up to T c is also of 
considerable interest. Especially important is to trace 
the smearing of the anomaly AC by the same effect of 
disorder which leads to the transformation of the intrin- 
sic exponential low-T behavior of C s (T) into a power-law 
one. These objectives were met by numerical calcula- 
tions. 

For this purpose, a Gaussian model distribution func- 
tion fo{So) was used: 



fo(S ) oc exp 



(*o ~ go) 
2d 2 



(12) 



The parameter <5q designates the peak position, which 
may vary from to I. By changing the parameter d 
we control the dispersion of the domain superconducting 
properties. Nevertheless, for any d the function /q (So) 
does not vanish in the limit Sq — and its Taylor se- 
ries begins with a constant as the main term. Only for 
highly improbable distribution functions, when simulta- 
neously fo(So) extends to So = and matches the con- 
dition fo{So = 0) = 0, the Taylor series may begin with 
the next term resulting in the asymptotics C S (T) ~ T 3 . 

In Fig. [l] the dependences (c(t)) are depicted in the 
panel (a) for Sq = 1 and different dispersion values d. 
A substantial spreading of the anomaly AC readily seen 
in Fig— [U-seems quite natural in view of the results for 
MgB2c2rt3. However, the concomitant superposition of 
various domain contributions distorts the whole curves 
C S (T) and C(T), which is much less trivial. This very 
superposition leads for low T to the power-law behav- 
ior, the asymptotics of which was analyzed above. The 
low-T parts of the curves (c(t)) are displayed on the log- 
log scale in the panel (b). Dotted straight lines corre- 
spond to the pertinent T 2 -asymptotics for each curve. It 
is clear that the validity range of the asymptotics extends 



with the increase of d. Although intervals where the T 2 - 
approximation holds good exist for any d, for small d it is 
merely of academic interest, because both temperatures 
and heat capacities become too tiny to be experimentally 
significant. On the other hand, for higher T in this case 
the averaged dependences (c(t)) lie rather close to the ex- 
ponential curve inherent to the BCS theory (the dashed 
curve). Such transitional parts of the dependences (c(t)) 
describe well the—exponential low-T behavior for some 
samples of MgB2L 3 Jtj with smaller exponents than in the 
BCS case. 

For large d, when the Gaussian distribution function 
/cP(<M becomes almost uniform (such a random dense, 
although quasi-discrete- distribution of gaps was found 
in point-contact spectraQ) , the quadratic asymptotics are 
valid at least up to t = 0.1 (for the uniform distribution 
fo (So) = const the relative error of the t 2 -asymptotics is 
« 0.6% at t = 0.1 and_£« 5% at t = 0.2), which agrees 
with the measurements^!!. For intermediate d the experi- 
mental data in the relevant T-range may be satisfactorily 
represented by power-law curves C(T) ~ T n with n > 2. 

One can make another important conclusion from the 
numerical data shown in Fig. [l]. A one-parameter fitting 
explains both the smearing of the heat-capacity anomaly 
at T c and the appearance of the power-law asymptotics. 
The latter reproduces the results appropriate to sttBer- 
conductors with order parameters of rfa^-HnVwavell^l or 



extended s-wave with uniaxial anisotrop; 



symmetry. 



The patterns displayed in these figures explain well the 
experimental heat capacity dependences C(T) for MgB2, 
which demonstrate power-law behavior for lowest_aitain- 
able THE! or above the exponential low-T taiO. At 
the same time, the reduction of the anomaly AC at T c 
with the increase of d, traced in Fig. |l|(a), adequately 
describes the AC magnitudes inferred from the analy- 
sis of the observed total heat capacity of MgB2, making 
allowance for crystal latt ice and papurity components. 
Namely, fi « 1.13EJ, 0.82BEI, 0.7E3, so that the exper- 
imental specific heat jump is substantially smaller than 
the BCS value pbcs- 

To summarize, we presented a phenomenological model 
of the disordered s-wave superconductor with a random 
domain network possessing continuously varying super- 
conducting properties. The characteristic domain size 
L is supposed to exceed the superconducting coherence 
length £o- The spatially averaged electronic heat capacity 
(C(T)) is calculated. It is shown that its low-T asymp- 
totics is a power-law one ~ T 2 , whereas the anomaly 
AC at T c is simultaneously smeared. These are just the 
features appropriate to the heat capacity of MgB2. Al- 
though we do not know the exact nature of the partition 
into domains in this compound, its very existence un- 
doubtedly manifests itself in Raman, point-contact and 
tunnel spectraB~H. One may speculate that the multiple- 
gap superconductivity originates from some kind of a 
phase separation rather than from the existence jof sev- 
eral groups of current carriers in the same volumeu. The 
presented theory may be also invoked to explain low- 
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T properties of cuprategla, although the microscopical 
background of the multi-gapness may be quite different 
in both cases. 
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FIG. 1. (a) Temperature dependences of normalized to- 
tal electronic heat capacity (c(i)} in comparison with the 
BCS-dependenceof superconducting phase fraction. Gaus- 
sian distributions with <5g = 1. (b) Low-temperature por- 
tions of the relevant curves on log-log scale together with their 
t 2 -asymptotics. 
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